perm filename FMSIG.F4[1,MUS]1 blob sn#088194 filedate 1974-02-19 generic text, type C, neo UTF8
COMMENT āŠ—   VALID 00002 PAGES
C REC  PAGE   DESCRIPTION
C00001 00001
C00002 00002		DIMENSION W1(512),FUNC(512),SU(100),ST(550)
C00004 ENDMK
CāŠ—;
	DIMENSION W1(512),FUNC(512),SU(100),ST(550)
	GO TO 60
40	CALL HYDPOG(2)
60	TYPE 50
50	FORMAT(' TYPE A,CF,MF,I1,I2,STEPS'/)
	ACCEPT 100,A,FC,FM,XI1,XI2,STEPS
100	FORMAT (6F)
	CALL TYPLOC(-300,-512)
	CALL DPYSET(1,SU,100)
	CALL DPYBRT(1)
	CALL AIVECT(0,0)
	CALL ALINE(0,-12,0,12)
	CALL ALINE(256,-12,256,12)
	CALL ALINE(-268,128,-244,128)
	CALL ALINE(-268,256,-244,256)
	CALL ALINE(-268,-128,-244,-128)
	CALL ALINE(-268,-254,-244,-254)
	CALL ALINE(-264,0,256,0)
	CALL ALINE(-256,-256,-256,256)
	CALL DPYOUT(1)
	XI=XI1
	XI1=(XI2-XI1)/STEPS
	N=STEPS+1
	CALL DPYSET(2,ST,550)
	CALL DPYBRT(2)
	Y=6.2832/512.
	DO 400 N1=1,N
	CALL CLRPOG(2)
	I=-1
	DO 200 J=1,512
	X=J-1
	FUNC(J)=A*(SIN((FC*X)*Y-XI*(COS((FM*X)*Y))+XI))
	IY2=FUNC(J)*256.
	IF(I.EQ.0)GO TO 150
	CALL AIVECT(-256,IY2)
	I=0
	GO TO 160
150	CALL SVECT(1,IY2-IY)
160	IY=IY2
200	CONTINUE
	CALL DPYOUT(2)
	XI=XI+XI1
400	CONTINUE
	TYPE 403
403	FORMAT(' TYPE 0 TO CONTINUE'/)
	ACCEPT 300,M
300	FORMAT (I)
	IF(I.EQ.0)GO TO 40
	END